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Abstract 

We present a complete and detailed derivation of the finite lattice 
spacing corrections to staggered fermion matrix elements. Expanding 
upon arguments of Sharpe, we explicitly implement the Symanzik im- 
provement program demonstrating the absence of order a terms in the 
Symanzik improved action. We propose a general program to improve 
fermion operators to remove 0{a) corrections from their matrix elements, 
and demonstrate this program for the examples of matrix elements of 
fermion bilinears and Bk- We find the former does have 0{a) corrections 
while the latter does not. 



1 Introduction 



With the ocuurence of a new generation of teraflop paraUel supercomputers, we 
wiU be able to simulate lattice QCD with smaller and smaller statistical errors. 
It is now of increased importance to gain control of various kinds of systematic 
errors which cither affect the numerical results directly or affect the way in which 
physical quantities are extracted. One of the most important systematic errors 
comes from the finite lattice spacing a which generates errors of the order of 
i^QCD- For the present lattice computation, this corresponds to the corrections 
of the order of 20% ~ 30%. Another important systematic error comes from 
the choice of lattice operators. There exist a variety of lattice operators which 
approach the same continuum operator in the limit a ^ 0. However, many 
of these operators differ from the continuum limit at 0(a), and a systematic 
formalism is needed to improve the lattice operators so as to remove these 0{a) 
corrections. 

For the case of Wilson fermions, the standard lattice action differs from 
the continuum quark action by a term of 0(a). So, both the action and the 
operators need corrections at order a. Applying the improvement program of 
Symanzik[0 to Wilson fermions, a procedure was proposed in ref.Q |^ to reduce 
the systematic errors due to the finiteness of the lattice spacing, from terms of 
0(a) to ones of 0{gQa), and it was numerically demonstrated in ref. Q| that this 
procedure can reduce the finite a corrections from 30% to around 5%. 

The meaning of the statement that there is no term of order a in the stag- 
gered fermion action is not clear. Let us use the free staggered fermion action 
as an example: 

Sf = ^a^x{x)rl^,{x)-^[x{x^ ^i) - x{x - ^ m'^ x{x)x{x) ■ (1) 

Following Golterman and Smit we denote the Fourier components of the 
fields X and x as x and Xi and decompose momentum space as 

k^p + TTA, (2) 

where fc^ S (0, 27r/a), e (0,7r/a), {TTA)fj. — A^-k, in which = 0, 1. If we 
define the fermion fields as 

'^(?') = 1E(-1)'''''^^*(p + ^b), (3) 



A,B 



where 



iA = it'itii'it\ (5) 
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we can write the action as: 

Sf = '^(P'>(^ "^A'" sinp^^a + m)tp{p). (6) 

p /J 

where fl is the lattice volume. It is clear that there is no order a term in 
the action, hence the free staggered action is accurate to 0{a^). However, the 

coordinate fields corresponding to ip (ip) are non-local superpositions of the x's 
(x's) over all the lattice sites. 

On the other hand, if we define the local hypercubic fermion fields as in 
ref.rtni 



A A 



where 



x = y + A, y^ = 0,±2,... (9) 
then the action in the momentum space can be written as: 

Sf = — ^^^(p){^[(7A««)/)^sinp^2a + a(75(g)^5;,)(-sinp^a)^]+m}g(p). 
p M 

(10) 

It is obvious that there are order a terms in the action and in the propagator 
for the fields q and q. In this case, we say that the fields q and q need to be 
improved. There exists a set of improved fields 

xi(y) = (l-a5]A^a^)XA(2/), (11) 

x'Aiy) = XA{yKi-aJ2A,Xl (12) 



where 



such that 



^f:fiy) = ^[fiy + 2^,)-f{y-2^i)], (13) 



= ^ E 9'(p)(E >^ ^i'^^''^" + "^)9~'(^^) + o(°')- (14) 

p p 

Note those improved fields are still local and superior to the nonlocal fields 
both computationally and theoretically when gauge couplings are included. So, 
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if we use the improved fields which remove the order a terms from the action to 
construct a lattice fermion operator, there will be no 0{a) corrections to its free 
field matrix elements. For Landau gauge, Sharpe ^ proposed the following 
smeared operator: 

Xa(j/)"^'' ^lY.^A{y + 2u[l - 2A,]). (15) 
It is easy to show that 

XA(y)™^'' = xi(y) + 0(«')- (16) 

The full staggered fermion action including gauge couplings is much more 
complicated. In section 2, we will give a set of improved fermion field variables 
in terms of which the action has no explicit order a terms at tree level. In section 
3, we will expand upon the argument given by Sharpe in ref Q to prove that 
there are no 0{a) terms which can be added to the staggered fermion action. 
Based on these two arguments, we conclude that staggered fermion action is 
already accurate to O(a^), and that we should use the improved field variables 
to construct fermion operators to reduce order a corrections from their matrix 
elements. We apply this program to the case of < 0|s754(i|i^° > and Bk as 
examples. We will also determine the additional operators that must be added 
to improve the standard staggered fermion currents to define operators whose 
matrix elements are accurate to O(a^). We list the lattice symmetry transfor- 
mation properties of the fermion fields needed in this paper in the Appendix 
A. 



2 Improving the fermion fields 

Working with the even position variable y of Eq. (^) , the lattice gauge-covariant 
derivative is defined as follows: 

Dfjiy) = j^[U^{y)U^{y + pi)f{y + 2pi)~Uliy~fi)Ul{y~2fi)f{y-2fi)]. (17) 

Following ref.||l^, we define the gauge covariant fermion fields as 

<fA{y) ^UA{y)xA(y), 

^A{y) = XA{y)u\{y), (18) 

where lAA{y) is the average of link products along the shortest paths from y to 
y + A. The classical continuum limit of action can be written as: 

SfVp, 0\°~^ I Y1 '?A{y){Y^ {it. ® I)AB^f^ + ^'^^AB ~ (75 ® ^5il)abDI 

•'y AB M p 
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+0(a'). (19) 

where = 9^ +117^^1 is the continuum covariant derivative. At first sight, the 
action contains order a terms. Likewise, it is clear that the fermion propagator 
for the hypercubic fields ip and (p deviates from the continuum propagator by 
terms of order a. However, if we introduce the following improved field variables 

xi(y) = (l-a^A,I?^')xA(y), (20) 
XA{y)=XA{y){l-aJ2A.Dl), (21) 

V 

and replace x, x in Eq.|l^ by and x^j then the classical continuum limit of 
the action can be written as 

-'y AB ^L 

We see the Feynman rules of the improved fields differ from those of the con- 
tinuum theory by terms of O(a^). In ref.Q it was proven that if the matrix 
elements of the tree level improved operators differ from the continuum ones 
only by terms which are of O(a^), then in the full theory there are no terms 
of 0{a) nor terms of O(a(7o"lii"a). We restate their argument here. Since the 
lowest order correction to the propagator and vertices are of order a^, a term 
proportional to a^Q^ln^a from an n-loop diagram can only occur if one of the 
loop integrals diverges like 1/a. If this is true, there will be at most n — 1 loga- 
rithmically divergent loop integrals which will result in a term of O(a(7o"ln"~^a) 
which behaves like 0((;§a) as a ^ 0. So we conclude that there are no terms 
of 0(a(7Q"'ln"a), and therefore the matrix elements differ from their continuum 
ones by terms at most of 0{gQa). 

Using the new fermion fields, we can construct improved fermion operators. 
For example, the improved fermion bilinears have the following form: 

XAiy)ilS «) ^F)ABXBiy) = XAiy){lS ® £.F)ABXB{y) 

-lfYdi[XA{v){{lS®S,F)AB - {lbvS®(,bvF)AB)XB{y)] 

V 

-7^^XA{y)\{l^vS ® ihvp) AB - {iSbu ® S,F5^) AB\DuXB{y) 

+0(a'). (23) 
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3 Improving the staggered fermion action 



In contrast to the calculation of matrix elements, the action is already accu- 
rate through order a to all orders in ^q, as we will now discuss. Thus physical 
quantities that depend only on the form of the action (for example, particle 
masses determined from correlation functions) will have no corrections of order 
(7o"a. This can be demonstrated by recognizing that if there were a correction 
of order ggU, we must necessarily be able to add some dimension-5 operators 
'^^i'^i9o^i^^ which must be invariant under the lattice symmetry transforma- 
tions to cancel this order ^qO correction(see ref. [Q Q Q). However, following 
ref.[^, we will now prove that there exists no dimension-5 operator (for the defi- 
nition of the dimension of lattice operators, see ref. [Q) which is invariant under 
the lattice symmetry group (rotations, axis reversal, translations, U{1) (E) U{1), 
charge conjugation), and therefore, no order a term can be added to the stag- 
gered fermion action. 

Following standard notation, we rewrite the staggered fermion action as 

y.y' A,B 

+mSiy - y')SAc]xc{y')''Uiy + A,y' + CT", (24) 

where ^ 

(is «> iF)AB = -^Trij^jsjBlp), (25) 

D^{y,y')AB ^D^{y,y')6AB + aAf,{y,y'){-f^r, (S) ^^5)^^, (26) 

in which 

D^.{y, y') = ^[5{y + 2^i- y') - 5{y ~ 2fi - y% 

A^{y, y') = ^ [S{y + 2fi ^ y') + S{y ~2^l~ y') - 25{y - y% 

For convenience, we will not write out the SU (3) links explicitly in the remainder 
of this section unless there would otherwise be confusion. Given an operator, 
the reader can write out the full form very easily. For example, starting with 
the operator X75 ® Cs^^^X would construct the corresponding gauge 
invariant operator by the following substitution: 

J2 J2 XAiy)h5 ® e5M)AB^M(2/. y')BcD^{y'y")cDXD{y") — > 
y',v" M 

E E XA{y){H^®I^)ABU{y + A,V + B)D^{y, y')BcU{y' + y" + C) 
y',y" A' 

xD^iy', y")cDU{y' + C, y" + D)xD{y"). (27) 

where U{y -I- v4, y + B) is the average of the products of link matrices corre- 
sponding to each of the shortest paths from point y + Aio y + B. 



5 



Using the transformation properties of the staggered fermion action (see Ap- 
pendex A in detail), we try to construct all symmetrical dimension-5 operators 
which have the general form xis ® £,Ff{D, D)x where / is a homogenous real 
polynimial of degree 2. 

Invariance under requires that S" + F is odd, so only the following 

combinations of (g) _F are valid: 

(/,75,7m-'75m.)® (6,^5a), (28) 

(7^,75/.)® (/,C5,CAr,C5Ar). (29) 

Under reflection with respect to a hyperplane normal to the p direction, we 
have the following transformation: 

X ^ IpX, (30) 

X~-Xl;\ (31) 

and 

D^^{l-26^p)IpD^I-\ (32) 

A^^XpA^X;!, (33) 

D^^{\~25^,p)IpD^I-\ (34) 

Using the transformation formulae of 75 ® listed in Eq.( ^9|) of the Ap- 
pendix, we deduce that axis reversal invariance and Ua{^) invariance allow only 
the following terms: 

75®^M^', (35) 

75®^5,.i?^, (36) 

lz[l^.,lA®^B{^^.+^•^)[D^.,D,], (37) 

l5[l^..lA'^U^^.~^.){Dp,D,}. (38) 

where Dp,D^ can be replaced by D^^^Di, without affecting these operators up 
to order . 

Under a rotation around the center of a hypercube, we have the following: 

X ^ n^'^h. (39) 
X ^jn^'"^~' (40) 
Dp ^ n'-p^^ Rp^D.n^p"^-^ (41) 

A^ -> 7^(''") A,7^(^")-l (42) 
Dp n'^P^^Rp^D.n'^P''^-^ (43) 

(44) 
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Combining the transformation properties listed in Eq. ( |8^ ) of the Appendix, 
the rotational invariance will further eliminate the term in Eq.(p5|) but allows 
the remaining three terms Eq.(|3^ - |3^ ). 

Finally, let's discuss invariances under translation by one lattice unit: 

X^S^P^X, (45) 
X^_xS^'^-\ (46) 
D,,^S^P^D,,S'-P^-\ (47) 

^ S^P^A^S^p'>-\ (48) 
D,^S(P^A^^%jD,S(P^-\ (49) 

where 

A^ilAy, v')ab = siF) f^SAB6{y - y') + a[{-l f^ (-1)^"] 

x[a5ABApiy,y') + J^i^^^^^^gDp{y,y')]'j , (50) 

4i~Fiy, y')AB = e{F) {{~ir^6AB6{y ~ y') + a[(-l)^'' - (-1)^"] 

x[a5AB'Ap{y,y') - [izp ® S.5p) ab^^ p{y,y')]\ , (51) 

and 

S^P^-^-is ® ^fS^"^ = js<»^FAiip. (52) 

From these properties, we can see that none of the terms listed in Eq.(^ - 
^8| ) are invariant under lattice translation! So, we conclude that there is no 
dimension-5 fermion operator which is invariant under the lattice symmetry 
group, and therefore no dimension-5 operator can be added to the staggered 
fermion action. 



4 Applications: 

As we argued above, actual numerical simulation should use the improved 
fermion field variables. However, in most situations, we can use the improve- 
ment program proposed in this paper to remove the 0{a) corrections without 
increasing the computational work. Here, we apply this program to the calcu- 
lation of the matrix element < 0|s754(i|iir° > which gives fx in the continuum, 
and the calculation of Bk ■ We will show that the former differs from its contin- 
uum counterpart by Olnika), but Bk has no 0{a) corrections. We also apply 
the improvement program to the matrix elements of lattice currents. 
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4.1 < 0|s754C^|i^° >: 

The axial current used in the (Landau gauge) numerical simulation is: 

My) = XAJn^i^^i^ABXBiv). (53) 

AB 

From the continuum expression 

P(t)™"* =< 0|A4(t)™"*|if" %/2//fmfee-'"''l*l, (54) 
we define, on the lattice, 

P{t)=<Q\Y^A^{x,t)\K° >, (55) 

and put the wall source that creates the on the time slice at t = 0. Then 
we will have 



where 

/± = /K±0(mfca). (57) 
If we don't consider 0{g'^a) terms, we can take only the term 



in Eq. ( |2^) because other terms contribute zero "flavor" trace at the tree level. 
So, we have 



^/aZ-^^^m.e-"^!*!, (t < 0) 



(58) 



and 



So, we get 



fV"'' = {l + \mka)f+ + 0{a^), (59) 
/K''"'' = (l-^m,a)/^ + 0(a^), (60) 



/I - (1 T ^™fca)/K + O(a'). (62) 
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10.0 15.0 20.0 25.0 

Time 



Figure 1: The value of ln(< 0|^4(t)|-ft'° >) with respect to the time t. Where 
ruda = irisa = 0.01. Calculated with the cubic wall source method. 

The numerical data (from the full QCD simulation on a 16'^ x 40 lattice with 
the cubic wall source, at the quark mass m^a = m^a = 0.01, see ref . |[T^ ) for the 
unimproved and improved matrix elements are shown in Figure |l|, from which 
we see that (/^ - « m^a ^ 25% and (/^'""^ - ^ 5% is 

much smaller. 

From this simple example, we can see that if we do not consider 0{g^a) 
corrections, the improved operator is eqivalent to the extrapolation: 

P{t) = V2/^"Pmfce-"'=l*+i/2|, (63) 

4.2 Bk: 

The formula for calculating Bk is: 



Mk 

where 



Bk-^. (64) 



Mk -< i^°|s7M(l + 75)ds7^(l + 75)d|^° >, (65) 
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K0\sj4rod\0 >< 0\s-f475d\K° > . (66) 
The improved numerator is (omitting the 0{gQa) terms): 

= Mk- \dlMK + O(a'). (67) 

Since M.K(f) is a plateau (i.e. time independent) within the statistical error, 
there is no 0{a) corrections to the numerator. 
Since 

/^/a- = /i + 0{a\ (68) 
the denominator also has no 0{a) corrections even if no attention is paid 
to an accurate definition of jx- Hence, we showed that there is neither 0(a) 
nor O(a(7o" log" a) corrections to Bk- Sharpe has examined this question in 
greater detail and argued that in fact there are no corrections of 0{g^a) also. 
However, if we calculated the denominator only in one time direction and took 
the square of (or /^), there would be an error of order of 0{mka). 

4.3 Renormalization of lattice currents: 

The lattice currents can be written as 

jLt = XlT^X, (69) 
and according to rcfp3[| , their renormalized continuum forms can be written as 

where Zj is the usual (divergent) renormaliation constant and Kj is a finite 
lattice renormalization constant. Using the method developed in this paper, we 
can explicitly determine the improved currents accurate to 0{a?). For exam- 
ple, the conserved vector current and axial vector current corresponding to the 
C/y(l) (g) 1/4(1) lattice symmetry can be written as follows: 

{y) = V,,{y) -^Y. dAxAh;:Wl)ABXB] 
a 



a \ ^ 



^df,[xAh^® I)abXb] 

,2\ 



^^d,[XA{l5M^^5.)ABXB] + 0{a'), (71) 



= 4'(y) - ^ E d.[XA{-/,5 ® i5)ABXB] 

a 
~2 



^^,[xA{l^I.b®ib)ABXB] 



'J E 9AxAh[^.M ® ^.)abXb] + O(a'). (72) 
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The effect of the second term on the right hand side is to shift the position 
y, labehng the current, from the corner to the center of the hypercube. The 
third term whose effect is to shift in the ^'s direction occurs here because the 
currents are non-local operators which involve an overlap between two nearest 
hypercubes. The forth term is a mixing of a different spin-flavor operator and 
is necessary to remove all order a effects from a general matrix element. 

5 Summary 

In this paper, based on the demonstration of that there is no dimension-5 
fermion operator which is invariant under all lattice symmetry transformations 
and that there exists a set of improved fermion fields with respect to which the 
tree-level action has no order a terms, we concluded that the staggered fermion 
action is already in fact improved to O(a^). We argued that to remove order 
a corrections from the matrix elements, the first step is to use the proposed 
improved fermion field variables to construct fermion operators so that they 
differ from the continuum ones by order of at the tree level. Then we showed 
that to all orders of perturbation theory treating gQ\na ^ 0(1), the correction 
is at most of 0{gQa). Furthermore, for the on-shell quantities, they are accu- 
rate to O(a^). We applied our program to the matrix element < 0\sj54d\K'~' > 
and found that the unimproved one differs from the continuum by a factor of 
0{mka). At the same time, we showed that there is no 0{a) corrections to Bk, 
which is consistent with the result of Sharpe Q . We also discussed the matrix 
elements of the lattice currents, and obtained the explicit terms which should be 
added to the original current operators to define improved operators accurate 
through 0{a). 



6 Acknowledgements 

I warmly thank Norman H. Christ who proposed this topic, stimulated me 
to finish this paper, and gave me the benifit of extensive discussions during 
every stage of this work. I am also grateful to Weonjong Lee for the numerical 
data appearing in this paper and many helpful discussions. I also thank Bob 
Mawhinney for the helpful and interesting discussions on the lattice symmetry 
of staggered fermions. Adrian Kahler and ChengZhong Sui checked out some 
typing mistakes, their kind hospitality is gratefully acknowledged. 

A Symmetry Properties of Staggered Fermion 

For completeness, we collect some formulae connected with the transforma- 
tion properties of the staggered fermion under the lattice symmetry group from 
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ref.Q. 

A.l U{1)a: 



XAiy) ^ e^'^^^hAiy) 

where 

e{A) = (-l)^*'^'' 
A. 2 Reflection with respect to a Hyperplane: 

The transformation of the fermion fields are: 

XAiy) ^^^ip{v,y')ABXB{v')-, 

B y' 

XA{y) ^^^XB{y')^p^{y',y)BA, 

B y' 

where 



ip(y, y')AB = (7p5 «) ibiAB^ihv - v')^ 
^p^iy,y')AB = {ibp®s.b)AB^i^py-y')^ 

and 

^ "^^^ " I -yp. (m = P) 

The spin-flavor matrices transform as: 



{is ® S.f) (7p5 C5)(7S ^f)(75p «> 6), 

which are explicitly listed as following: 



5* 








I 


I 




I 


75 - 


* -75 


?5- 


>C5 


7p - 


. (1 - 25^p)i^ 


^p- 


>-^P 


75m - 


-(1 - 25^p)75^ 


C5m ^ 






(l-2(5^p)(l-2J,p)7^, 


Cpi' 




lhp,u 


^-(l-2,5pp)(l-2^,p)75^. 




^ ^hp,v 
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A. 3 Rotations by 7r/2 arround the center of a hyperplane: 



Define R^ii, such that 



Vcr, (M = P) 

{Ry)ix = RiJ^uVv = { -Vp, (m = 0-) 

Then, we have the following transformation: 

XA{y) ^^Y.n^'"'\y,y')ABXB{y'), (80) 

B y' 

XA{y) --Y.T.^Biy')Ti^''^Hy,y')A'B, (si) 

B V' 

where 

n'-P''\y,y')AB = liil -lp.)®{^a-Cp))AB^iR''y - y'), (82) 

n^'"^ {y, yTAl = \W+i^^)®(^^^)AB5{Ry - y'). (83) 

The spin- flavor matrices transform as: 
1- 



[IS ® ^f) ^{{l - Ipa) ® {^a - ip)){lS ® ^f) 2 ((1 + Tpa) ® {^a " ip))- (84) 

which are explicitly listed as following: 

S F 
I ^ I I 

iij.-^Rp\i\ ^p^-\Rp\\£,x 

75At RpX75\ S.5p — > \Rpx\(5X 

Ipv RpxRurlXr ^pu ^ \Rpx\\RiyT\i\T 

Ibpv — > RpxRvrlbXr ^5pv — > —\Rpx\\Rut\^5Xt 

A. 4 Translations by one lattice unit: 



Tp-- 



Xp — Xp + 1, 

x'p = Xp, in i- p) 
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lead to the following transformation on the XA, XA^s: 

XA(y)-EE'5*'"(y,y')ABXB(y'), (85) 

B y' 

xA{y) - EE^^(y')5(''Hy,2/')AB, (86) 



where 



B y' 



S^P^ {y, y')AB = \ {{TWlp - 7p5®^5)AB^(y - v') 

+ (T¥T; + 7^;r¥I^)Ai3<5(y + 2p-y')], (87) 
5^") (y, v')ab = \{(j®ip ~ iW®1^)abKv - 2p - y') 

+(/ ®ip + 7p5 ® i^)ABKy ~ y'% (88) 
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